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Abstract 

The temperature dependence of the giant magnetoresistance (GMR) for currents par- 
allel and perpendicular to the multilayer plane, is discussed by taking account of the 
random exchange potentials, phonon scatterings and spin fluctuations. The effect of spin 
fluctuations, which plays an important role at finite temperatures, is included by means 
of the static functional-integral method developed previously by the present author. Our 
model calculations well explain the observed features of the parallel and perpendicular 
GMR of Fe/Cr and Co/Cu multilayers recently reported by Gijs et al. 
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I. INTRODUCTION 

The giant magnetoresistance (GMR) 1 in magnetic multilayers is one of the most at- 
tractive subject in current solid-state physics. In recent years much progress has been 
made in understanding the GMR and its related phenomena. 1 One of the important as- 
pects of the GMR is its temperature dependence. A careful study of the temperature 
dependence of GMR is not only important in understanding its mechanism but also very 
useful to its realistic applications. Most of the magnetic multilayers are fabricated with 
transition metals such as Fe, Ni and Co. It would be instructive to briefly discuss the 
temperature dependence of the resistivity of bulk transition metals, before we study the 
temperature dependence of the GMR or of the resistivity of transition-metal multilayers. 
It has been reported that when the temperature is raised from T = K, the resistivity 
of Fe 2 or Co 3 gradually increases up to the Curie temperature, where it has a cusp (see 
Fig.l). This characteristic temperature dependence of the resistivity is interpreted as due 
to the contributions from impurity, phonon and magnetic terms. The last contribution is 
classically discussed as spin-disorder scatterings with the use of the s-d model. 4 Lately, a 
modern theory on the itinerant-electron magnetism has accounted for it in terms of spin 
fluctuations. 5 

It has been well known that d-electrons in transition metals show both the localized 
and itinerant character: the Curie- Weiss susceptibility and the large specific heat peak 
near the Curie temperature are easily explained by the localized-spin model whereas the 
non-integral ground-state moment and the large linear-specific heat coefficient favor the 
band model. It has been realized that the effect of spin fluctuations plays essential roles 
to reconcile the duality of d electrons. 6 The finite-temperature band theory, which has 



3 

been proposed by Hasegawa, 7 includes the effect of spin fluctuations by means of the 
static functional-integral method combined with the coherent potential approximation 
(CPA). Spin fluctuations including spin waves are shown to yield the T 2 contribution 
to the resistivity at T ~ by several approaches. 8 This type of theories 8 is, however, 
valid only at very low temperatures. In our finite-temperature theory, 7 spin fluctuations 
are regarded as localized, static modes with the adopted approximations. This method 
has proved useful in understanding the overall finite-temperature properties of transition 
metals, alloy and multilayers, 9 covering both below and above the Curie temperature. 

By employing the finite-temperature band theory, 7 we discussed in previous papers 10 
the temperature dependence of the MR ratio for currents parallel to the multilayer plane. 
The observed temperature dependences of Fe/Cr, 1M,n NiCo/Cu, NiFe/Cu and CoFe/Cu 12 
multilayers have been shown to be well explained by our theory. It has been pointed out 9c 
that a multilayer in which the normal and inverse GMR 13 coexist, may have an interesting 
temperature dependence beneficial for real applications. 

One of the purposes of the present paper is to generalize our theory 10 to the per- 
pendicular GMR, whose experimental 14-16 and theoretical study 17-20 has been currently 
performed. The other purpose is to include the phonon contribution to the conductivity 
calculation, which was neglected in our previous study. 10 The paper is organized as follows: 
In the Sec. II, we present our formulation applying our finite-temperature band theory to 
the GMR. Numerical calculations of the parallel and perpendicular GMR of Fe/Cr and 
Co/Cu multilayers are reported in Sec. III. Supplementary discussions are given in Sec. IV. 
II. CALCULATION METHOD 
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A. An Adopted Model and the Expression of GMR 

We adopt an A/B multilayer consisting of magnetic A and nonmagnetic B atoms with 
the simple-cubic (001) interface. The layer parallel to the interface is assigned by the 
index n (— 1 — Nf). The thickness of the A and B layers is assumed to be thinner than 
the mean free path and sufficiently thin compared with the spin diffusion length. 17 It is 
assumed that atoms A and B are randomly distributed on layer n with the concentrations 
of x n and y n , respectively (x n + y n = 1). The film is described by the single-band Hubbard 
model, in which the atomic potential (the on-site interaction) is assumed to be given by e A 
and e B (U A and U B ) when a given lattice site is occupied by A and B atoms, respectively. 

In order to study the finite-temperature properties of the magnetic film, we apply the 
functional-integral method within the static approximation to the Hubbard Hamiltonian. 7 
The partition function is evaluated by calculating the partition function of the effective 
one-electron Hamiltonian including the random charge and exchange fields with the Gaus- 
sian weight. The charge field is include by the saddle-point approximation and the ex- 
change field by the alloy-analogy approximation with the CPA. The energy-dependent 
coherent potential for an s-spin electron (s =t, I) on the layer n, Yi ns (e), is determined 
by the CPA condition. The coherent potentials, the average of the magnetic moments 
on the layer n, (M n ), and its root-mean-square (RMS) value, ((M n ) 2 ) 1//2 , are calculated 
self-consistently, details having been given in Ref.7. 

When we employ the CPA, the conductivity of the film is given by 20 

a(v = (I) (~rr) J d£ Tr (^ ImG ^ ImG ) (t,V = x,y,z), (1) 

provided the vertex correction is neglected. In Eq.(l) v% is the velocity operator and G is 
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the Green function matrix. The conductivities for currents parallel (||) and perpendicular 
(_L) to the film layer are given by 10,18,20 



with 




which is valid within the Born approximation. In Eqs. (2)-(5) A ns = Re E ns (e), A ns = |Im 
S ns (e) |, S ns is the coherent potential of an s-spin electron on layer n, and a n i s and u x 
are specified by the electronic structure of the film (see Eqs. (19) and (20) in Ref. [10a]). 
Analytic expressions given by Eqs.(2)-(5) have clear physical meaning. When currents flow 
parallel to the plane, an s-spin electron propagating successively from a site on layer n to 
a site on layer m, is scattered with the strength proportional to A ns and A ms , respectively, 
and its conductivity is given as a sum of such processes with the weight of a nms r nms . w 
On the contrary, in the case of the perpendicular current, the s-spin conductivity is given 
as of a series circuit of resistivities on successive layers, each of which is proportional to 
A„ s . 18 ' 21 In both cases, the total conductivity is a sum of the up- and down-spin channels. 
The so-called spin-flop process is implicitly included through the spin-fluctuation term 
which is responsible to a decrease in layer magnetization, as will be shown shortly. In the 
next section, we will employ our formalism in a semi-phenomenological way to discuss the 
temperature dependence of the MR ratio. 
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B. A Semi-phenomenological Study of GMR 

We adopt a system consisting of magnetic (Mi, M 2 ) and nonmagnetic (Ni, N 2 ) layers, 
whose thickness are M and N, respectively. Bulk scatterings are assumed to be important 
in these layers, related discussion will be given in Sec. IV. When magnetic moments on 
Mi and M 2 layers are in the antiferromagnetic (AF) configuration, the real and imaginary 
parts of the coherent potentials are given by 10 

A£f-iA£J = A s -iA s forneMi, (6) 

= A_ s - iA_ s for n G M 2 , (7) 

= A -«A forneNi,N 2 . (8) 

Using Eqs.(2)-(8), we get the parallel and perpendicular conductivities given by 10 

AF|| V"^ / 2c MM , C NN , a AF ( 1 , 1 \ , jAF f _}■ , \ \ , 

(9) 

^ = (i) 2 ^i:{ MAt+M l_ j+2jVAo }. no) 

with 

4m = Nf\e/h) 2 J E a nmT nm , (11) 

ngMi mGM 2 

< F = Nj\e/hfJ ]T E a nmTnm, (12) 
n<EMi mgMi 

and c^n; c mn; an d ^n F are gi ven by similar expressions. We employed the T = limit of 
Eqs. (2) and (3) because the relevant temperature is much less than the Fermi energy, e F . 
In Eqs.(9)-(12) u x = u x (e F ), Nf = 2(M + N), and the spin dependence in a nms and T nms 
is neglected. Subscripts, MM, NN and MN, denote the contributions from the interlayer 
scatterings between magnetic layers, between nonmagnetic layers, and between magnetic 



and nonmagnetic layers, respectively. On the contrary, the single subscript, M (N), ex- 
presses the contribution from the intralayer scatterings within magnetic (nonmagnetic) 
layers. 

On the contrary, when magnetic moments on the subsequent magnetic layers are in the 
ferromagnetic (F) configuration, the real and imaginary parts of the coherent potentials 
are given by 10 



AL - iAl 



A, - iA, 



Ao - ?A C 



for n G Mi,M 2 , 



for n E Ni,N 2 . 



(13) 
(14) 



We get the parallel and perpendicular conductivities given by 10 

8c F 



' F F 
C MM _|_ C NN _|_ 



H F H F ' 
^ ^M _|_ _zi 



V l A s A (A s + A ) A. 

*Ms) Vw '?{2matW;}' 

The MR ratio, AR/R, is given from Eqs. (9), (10), (15) and (16), by 
/ AR\ X _ (R AF - R F ) _ (a F - a AF ) _ (a - (3) 2 vX 



V R 



R F 



rAF 



4af3 



with 



a = A T /A , p = A l /A , 



X 11 = 



1 1 

+ 



(15) 
(16) 



X x (A=||,±), (17) 



x- 



) + {m) U 



a + 1 ' P + 11 +92 {m) {a + ^ 



N\ 2 ( 1 N ' 



(18) 

i 

i 

(19) 
(20) 
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In Eq. (19) g , g\ and g 2 are defined by 



d>M 2cmn ( N 
= 9o, = 9i 



4c 



MM 



CMM 



M 



Cnn + djq 



2c 



MM 



N\ 

92 [m) 



(21) 



The expression for the GMR given by Eqs.(17) and (20) is just the same as that derived by 
Edwards et al. 21 using the resistor network model and has been employed for an analysis 
of (AR/R) ±U 

Setting iV = in Eqs. (17)-(20), we get 

'A*y_ (a-py (22) 

(23) 



R J 4a/3[l + g (a + [3) 2 /apY 



V R J Aafi ' 

and the ratio of the parallel GMR to the perpendicular one is given by 

(Ai?/jR) " = / 1 , 9/ = < 1 (a = a/P). (24) 

(AR/R) [I + 9o{a + 1) 2 /a] v /HJ K J 

When electrons flow perpendicular to the layer plane, all electrons pass through the 
adjacent two magnetic layers. On the contrary, it is not the case for currents parallel to 
the plane; some electrons go through only the one of magnetic layers without probing 
the other magnetic layer. The second g term of the denominator of Eq.(22) denotes this 
contribution. Itoh et al. 20 claim that the anisotropy of the velocity operator : v^jv^ > 1 
is the main mechanism leading to (AR/R)^ < (AR/R)- 1 . The factor, i/H or z/^, is not, 
however, relevant because it is cancelled out when the MR ratio given by Eq.(17) is 
calculated. 

The temperature dependence of the GMR arises from those of a and f3, which is 
expressed in terms of the coherent potential of the film (Eq.(18)), whose imaginary part 
in the magnetic (Mi or M2) layer is given within the Born approximation by 10 

A s = ^ + A: + A p s , (25) 
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with 




(26) 



(27) 




p m Ps z(T/e m ) 



(28) 



where e A and e B are the spin-independent Hartree-Fock potentials and p s is the density of 
states of an s-spin electron at the Fermi level. The first term (A r s ) in Eq. (25) arises from 
the scattering due to random Hartree-Fock potentials for an s-spin electron: the second 
term (A|) comes from the effect of spin fluctuations: the third term (Af) is introduced 
for phonon scatterings whose explicit form will be given shortly (Eq.(45)). 

On the other hand, the imaginary part of the coherent potential in the nonmagnetic 
(Ni or N 2 ) layer, is given by 



where the first and second terms denote the contributions from random potentials and 
phonons, respectively, and po is the density of states at the Fermi level of the nonmagnetic 
metal. In Eqs.(28) and (29) 9 m and O n are Debye temperatures, and P m and P n are 
related with the electron-phonon interactions in magnetic and nonmagnetic metals. 
Using Eqs. (18), (26)-(29), we get a and (3 given by 



A = A r + P n p Z(T/Q n ), 



(29) 



P = A 



a = A 



(1 + 7 (T)) [xy(B + m{T)f + x(fi(T) 2 - m(Tf) + p m Z(T/Q m )} 

(i+ Po z(T/e n )) 

(1 - 7 (T)) [xy(B - m{T)f + x(p(T) 2 - m(T) 2 ) + p m Z(T/Q m )] 

(l+ Po Z(T/e n )) 



(30) 



(31) 
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with 

m(T) = (M A )/M , (32) 

p(T) = y/((M*)*)/M , (33) 

7 (T) = (p T -p i )/(p T +p i ), (34) 

A = np(U A M /2) 2 /A r , (35) 

B = (2/U A M )(e B -e A ), (36) 

p m = P m /n(U A M /2) 2 , (37) 

p = P n po/AS = p m A (P„/P m ) (p /p), (38) 

where p = (1/2) (pf + pj_) and M is the ground-state magnetic moment. 

At T = OK where m(0) = p(0) = 1 and 7(0) = 70, Eqs. (30) and (31) become 

a = a (T = 0) =xyA{l + l0 ){B + l)\ (39) 

A) = P(T = 0) = xy A (1 - 7o )(S - l) 2 , (40) 
from which the coefficients A and B are expressed in terms of a , f3 and 70 as 



1 ( I a ■ / fa \ (41) 



4a^/ I V 1 + 7o V 1 - 7o/ 



The normalized magnetic moment, m(T), and its RMS value, p(T), are in principle 
calculated with the use of the finite-temperature band theory. 9 We here, however, adopt 
simple, analytic expressions of m(T) and fJ>(T) for our model calculation, given by 10 

m(T) = y/l- (T/T c ) 2 , p(T) = l. (43) 
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The temperature dependence of the spin asymmetry j(T) defined by Eq.(34) is assumed 
to be given by 

7 (T)=7o m(T). (44) 

As for the phonon contribution given by Z(T/Q m ) in Eqs.(28) and (29), we adopt the 
simple Griineisen function: 

Z(T/6 m ) = (T/9 m ) 5 l &m/T dy l —^ r — y (45) 

which is 124.43 (T/6 m f at T/0 < 1 and T/46 m at T/0 m > 1. 

Now we may calculate the MR ratio, AR/R, as a function of temperature with the use 
of Eqs. (17), (19), (20), (30), (31), (41)-(45), when we treat a , (3 , 7o , g Q , g u g 2 , T c , G m , 
©n, Pm, Po, and y, as input parameters. Our strategy for calculating the temperature- 
and layer-thickness-dependent MR ratio is as follows: We first determine the parameters, 
a , Po and 7 to be consistent with the band calculation, and also g , g 1 and g 2 so as 
to reproduce the iV dependence of the observed, ground-state parallel GMR. Then fixing 
there six parameters thus determined, we calculate the finite-temperature GMR with the 
additional parameters, Tc, © m , n , Pm, Pn and y, which can be properly chosen, as will 
be discussed in the model calculations of the next section. 

III. MODEL CALCULATIONS 

A. Fe/Cr Multilyers 

Gijs et al. 15 have observed both the parallel and perpendicular GMR for a sample of 
(3 nm Fe + 1.0 nm Cr) multilayer, whose results are plotted by circles and squares in 
Fig. 2, respectively. 
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Firstly we consider the case of T = 4.2 K. We determine the value of 70 = 0.4 from the 
ground-state band calculation of Pt/Pl = 2. 3. 22 We adopt a = 7.9 and f3 = 1.0, leading 
to B — 3.38 (Eq.(42)), which is consistent with the value estimated from Eq.(36) by using 
the band parameters such as e Fc etc. We choose the parameters of g = 0.045, g\ = 0.77, 
and g 2 = 3.05, such that we have a good fit to the envelope of the observed layer-thicknes 
(t N ) dependence of parallel GMR in (3 nm Fe + t N Cr) multilayers. 11 

Next we consider the MR ratio at finite temperatures. We assume the Curie tem- 
perature of the multilayer of T c = 1000 K because the thickness of the Fe layers of the 
adopted Fe/Cr multilayers 11,15 is sufficiently thick to sustain the Curie temperature of 
bulk Fe. The Debye temperatures of Fe and Cr are assumed to be 6 m = 6 n = 460 K. 
The phonon parameters, p m and p , can be determined as follows: The total resistivity, 
R, of a pure, bulk metal is given from Eqs.(2), (3) and (30), by 

R(T) oc |^[(1 + s 70 m(T))(MT) 2 " m(T) 2 + p m ^(T/0™))]- 1 j , (46) 

from which the ratio of the phonon contribution, R p , to the total resistivity at T = T c is 
given by r p = R P (T C )/R(T C ) = p m Z(T c /& m )/[l+p m Z(T c /e m )}. The value of p m = 0.69 
is chosen from the experimental data of r p = 0.27 of bulk Fe (Fig. 1(a)). 2 We calculate p by 
Po — Pm ^ (Po/p) derived from Eq.(38) with P n = P m and p /p = O.7. 22 The parameters 
discussed above are summarized in Table 1. The solid curve in Fig. 1(a) expresses the 
resistivity, R(T), of bulk Fe calculated by using Eqs.(43) and (46) with 7 = 0.4 and 
p m = 0.69, which well reproduces the observed data. 2 

The last parameter y, which expresses a concentration of nonmagnetic atoms in the 
magnetic layer and which depends on a sample employed in a experiment, is treated as 
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an adjustable parameter. The parallel and perpendicular GMR of the Fe/Cr multilayer 
calculated with y = 0.002, 0.005 and 0.01 are shown in Fig.2. Our calculation with 
y = 0.005 well explains both the (AR/R)^ and (AR/R) L observed by Gijs et al. 15 

In order to study the temperature dependence of the GMR in more detail, we show in 
Fig. 3, A s (s =|, |) as a function of the temperature. When the temperature is raised, Af 
and Aj increase because of the contributions from spin fluctuations and phonons. Then 
the ratio, A 1 /A l (= a/ (3), changes from 7.9 at T = to unity at T > T c . Fig.3 also 
shows the decomposition of A s to various contributions from random potentials (A r s ), spin 
fluctuations (A*) and phonons (A§). We note that at T = T c , A 8 / A = 0.70, A p /A = 0.26 
and A s /A p = 2.65. This shows a significant spin-fluctuation contribution, as suggested 
from the resistivity data of bulk Fe. 2 

B. Co/Cu Multilayers 

We have performed a similar calculation to explain the temperature dependence of 
parallel and perpendicular GMR of the (1.2 nm Co + 1.1 nm Cu) multilayer observed by 
Gijs et al. ie We adopt a = 0.7, f3 = 8.4 ((3 /a = 14), 16 and 70 = —0.7 which comes from 
the ground-state band calculation of Pi/p[ ~ 0.15 of bulk Co. 22 We cannot determine the 
values of g , g\ and g 2 because the layer-thickness dependence of the parallel GMR of this 
series of samples has not been reported. Then we tentatively adopt g = 0.13, g\ = 0.39 
and g 2 = 0.11 by scaling the data of similar Co/Cu multilayer 23 as to reproduce the 
observed ground-state value of (AR/R)^ = 0.43. 16 The Curie and Debye temperatures 
are taken to be T c = 1400 K and Q rn = Q n = 445 K. We adopt p m = 1.62 from 
the observed ratio of r p = 0.56 for bulk Co (Fig. 1(b)), 3 and p /p = 0.3. 22 Adopted 
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parameters are shown in Table 1. The solid curve in Fig. 1(b) denotes the temperature- 
dependent resistivity of bulk Co calculated by using Eqs.(43) and (46) with 70 = —0.7 
and p m = 1.62. 

The calculated (AR/R)^ and (AR/R) 1 - of the Co/Cu multilayer are shown in Fig.4, 
where y is treated as an adjustable parameter. Both the parallel and perpendicular GMR 
observed by Gijs et al. 1& are fairly well explained by our calculation with y = 0.005. 

Figure 5 expresses the temperature dependence of A s and its components, A r s , A s s 
and Ap, which shows that at T = T c , A 3 / A = 0.42, A p /A = 0.53 and A s /A p = 0.79. 
Comparing these figures with the corresponding ones of Fe/Cr systems, we note that spin- 
fluctuation contribution in Co/Cu multilayer is less significant than in Fe/Cr multilayer. 
This fact is expected to be the main reason why the observed temperature dependence of 
the GMR in Co/Cu multilayer is less considerable than that in Fe/Cr multilayers. 

IV. CONCLUSION AND DISCUSSION 

We have discussed the temperature dependence of the GMR for currents parallel and 
perpendicular to the multilayer plane. We have included contributions from the random 
exchange potentials, spin fluctuations and phonons, which are considered to be main 
scattering mechanisms yielding the resistivity in transition-metal multilayers. Our model 
calculations have accounted for the following features of the observed GMR: 11 ' 15 ' 16 (1) 
both the parallel and perpendicular GMR are significantly temperature dependent than 
the (average) layer moment, (2) (AR/R) 1 - is larger than (AR/R)^, (3) the temperature 
dependence of (AR/R) 1 - is more significant than that of (AR/R)\ and (4) the temper- 
ature dependence of GMR in Co/Cu multilayers is less considerable than that in Fe/Cr 
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multilayers. The effect of spin fluctuations plays an important role to account for these 
three items whereas phonons play a secondary role. In fact, the items (l)-(3) can be 
explained without invoking phonons. 10 

In our phenomenological analysis, we have assumed that the bulk scattering is pre- 
dominant. On the contrary, when we take into account only the interface scattering, the 
expression for the GMR is given again by Eqs.(17)-(20) but with M replaced by /, the 
thickness of the interface, and with a and (5 expressed in terms of the quantities relevant 
to the interface. Then, they have ostensibly similar T and N dependence to those in 
which only the bulk scattering is included. It is possible to extend our analysis taking 
into account both the interface and bulk scatterings, although the calculation becomes 
laborious because it inevitably needs much number of parameters. Among many param- 
eters, the most important ones are go, oto, (3q and y; (AR/R)W generally becomes smaller 
than (AR/R) 1 - by go, and the essential feature of the temperature dependence of the 
GMR is determined by the ratio of ao/Po and y. 10,12 
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parameter 


Fe/Cr 


Co/Cu 


a 


7.9 


0.7 


A) 


1.0 


8.4 


7o 


0.4 


-0.7 


Po/p 


0.7 


0.3 


T c (K) 


1000 


1400 


e m , n (K) 


460 


445 


Pm 


0.69 


1.62 


go 


0.045 


0.13 


9i 


0.77 


0.39 


92 


3.05 


0.11 



Table 1 Parameters adopted for numerical calculaitons (see text). 
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Figure Captions 

Fig. 1 The temperature dependence of the observed resistivity (circles) of (a) bulk Fe 
(Ref.3) and (b) Co (Ref.4); the calculated resistivity, R(T), and its phonon term, R p (T), 
above Tc are shown by solid and dotted curves, respectively, results being normalized by 
R c = R(T C ). 

Fig. 2 The temperature dependence of the parallel (||) and perpendicular AR/ R (_L) 
of (3 nm Fe + 1.0 nm Cr) multilayers. Dotted, solid and dashed curves denote the 
calculated results with y = 0.002, 0.005 and 0.01, respectively; circles (squares) expressing 
the observed parallel (perpendicular) GMR (Ref.15). 

Fig. 3 The temperature dependence of A s of up-spin (solid curves) and down-spin 
electrons (dashed curves) calculated with y = 0.005 for the (3 nm Fe + 1.0 nm Cr) 
multilayer. Also shown are their decomposition to various contributions from the random 
exchange potentials (A£), spin fluctuations (A*) and phonons (A§); the calculated results 
being normalized by A c = A S (T C ). 

Fig. 4 The temperature dependence of the parallel (||) and perpendicular AR/R (_L) 
of (1.2 nm Co + 1.1 nm Cu) multilayers. Dotted, solid and dashed curves denote the 
calculated results with y = 0.002, 0.005 and 0.01, respectively; circles (squares) expressing 
the observed parallel (perpendicular) GMR (Ref.16). 

Fig. 5 The temperature dependence of A s of up-spin (solid curves) and down-spin 
electrons (dashed curves) calculated with y = 0.005 for the (1.2 nm Co + 1.1 nm Cu) 
multilayer. See a caption of Fig. 3. 



